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摘 要 : 气象 站 和 卫星 降雨 资料 估算 降雨 侵蚀 力 时 存在 无 法 反映 空间 异 质 性 且 精 度 差 的 问题 ,基于 
空 尺 度 对 中 国 的 降雨 侵蚀 力 进 行 评 估 , 并 结合 降雨 


CLDAS 多 源 融 合 降水 ,利用 lwo 模型 从 不 同 的 时 


量 \ 侵 人 包 性 降 十 次数、 侵蚀 密 度 等 指标 ,探讨 降 R 
雨 侵蚀 力 与 地 


= 


对 土壤 侵蚀 的 潜在 作用 。 结 果 表 明 :(1) CLDAS 降 
实测 数据 在 不 同 的 时 间 尺 度 均 有 良好 的 回归 关系 ,相关 系数 达到 0.8 以 上 ,与 


CMORPH 降雨 侵蚀 力 相 比 ,其 相对 误差 显著 降低 ,可 以 准确 反映 全 国 范围 的 降雨 侵蚀 力 季节 性 变 
异 。(2) 在 2001 一 2020 年 ,不 同 雨量 区 的 降雨 侵蚀 力 、 降 雨量 和 侵蚀 性 降雨 次 数 的 变化 趋势 基本 一 
致 ,高 雨量 区 的 年 际 变化 波动 剧烈 ,侵蚀 性 降雨 次 数 和 暴雨 过 程 协同 影响 降雨 侵蚀 力 的 大 小 。(3) 
空间 上 ,中 国 的 降雨 侵蚀 力 值 的 特点 为 东南 沿海 地 区 高 西北 内 陆地 区 低 。 时 间 上 ,侵蚀 性 降雨 集 


多 ,降雨 侵蚀 密度 越 大 。 
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文章 编号 : 


众所周知 ,降水 是 导致 水 土 流失 的 重要 原因 之 
一 ,其 主要 通过 雨滴 落 入 地 表 , 产 生 汇 流 冲 刷 表层 
土壤 ,从 而 对 地 形 和 土壤 结构 产生 破坏 ,最 终 导 致 
土壤 颗粒 物 和 有 机 物 等 营养 物质 随 雨 水 流失 "。 降 
雨 侵 刨 力 是 水 土 流失 重要 的 外 在 驱动 因素 ,可 以 反 
映 降 雨 强度 、 降 雨 历时 、 降 雨 频次 等 指标 对 土壤 侵 
蚀 的 贡献 程度 ,准确 估算 该 参数 能 够 为 区 域 土壤 侵 
蚀 危 险 程 度 的 评估 葛 定 基础 ,可 有 效 提高 对 区 域 土 
二 侵蚀 程度 预报 的 准确 性 。 

降雨 是 一 种 复杂 的 自然 现象 ,这 使 得 降雨 侵蚀 
力 的 直接 观测 成 为 一 大 难点 ,目前 对 于 降雨 侵蚀 力 
的 估算 主要 通过 地 面 气象 站 和 卫星 降雨 资料 结合 
不 同 的 模型 获得 ,而 对 于 降雨 侵蚀 力 的 分 析 则 主要 
侧重 于 在 全 国 尺度 或 局 地 尺度 进行 时 空 研 究 .趋势 
分 析 等 。 胡 琳 等 “利用 陕西 省 的 气象 站 点 资料 采用 
反 距 离 权 重 空间 插值 方法 绘制 了 全 省 的 侵蚀 力 空 
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中 在 5 一 8 月 , 夏 、 秋 两 季 对 土壤 造成 的 侵蚀 影响 更 大 。(4) 通过 对 年 降雨 量 、 年 侵蚀 密度 和 年 暴雨 
量 进行 分 区 定量 分 析 , 结 果 表 明 暴 雨量 与 侵蚀 密度 成 正 相 关 关 系 , 即 年 降雨 量 一 定 ,暴雨 事件 越 


间 分 布 图 ,结果 表明 降雨 侵蚀 力 在 不 同 季节 间 的 差 
异 较 大 ,陕西 北部 和 关中 地 区 主要 在 7 一 9 月 发 生 水 
土 流失 。 扫 水 清 等 “基于 气象 站 逐日 降水 数据 估算 
得 到 全 国 的 降雨 侵蚀 力 , 并 利用 区 里 金 插值 方法 得 
到 时 空 连续 的 结果 ,其 中 侵蚀 力 最 高 值 在 广东 东 
兴 , 为 2x10: MJ -mmhmmh -a FRED KEK AY 
蚀 力 最 低 , 仅 为 3 MJ mm:hm ”ha 。 伊 力 哈 木 5 
利用 38 个 气象 站 资料 分 析 了 新 疆 维吾尔 自治 区 的 
降雨 侵蚀 力 空间 变化 特征 ,降雨 侵蚀 力 总 体 呈 上 升 
趋势 ,与 降雨 量 的 空间 分 布 基本 一 臻 ,大致 呈 现 中 
部 高 .南北 低 的 态势 ,气候 和 海拔 高 度 是 影响 新 疆 
地 区 侵蚀 力 格局 的 主导 因素 。 

随 着 遥感 技术 的 发 展 , 卫 星 降 雨 产 品 为 降雨 侵 
人 蚀 力 的 计算 提供 了 新 的 思路 “。Vrieling 等 "利用 
TRAMM (Tropical rainfall measuring mission ) 降雨 数 
据 分 析 了 非洲 降雨 侵蚀 力 的 时 空 变 化 ,并 对 IMERC 
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(Integrated multi-satellite retrievals for GPM ) .CMORPH 
等 30 min 分 辨 率 的 卫星 降水 产品 在 降雨 侵蚀 力 方 
面 的 研究 作出 展望 。 王 凯 等 “使 用 订正 后 的 
TRMM3B42 降水 数据 计算 中 国 的 降雨 侵蚀 力 , 评 佑 
指标 表明 TRMM 侵蚀 力 在 不 同 的 空间 尺度 上 均 有 
良好 的 适用 性 。 黄 进 平 " 采 用 章 文 波 公式 "将 
IMERG 数据 应 用 在 甘肃 地 区 的 降雨 侵 刨 力 研 究 中 ， 
并 结合 随机 森林 进行 降 尺 度 ,验证 了 该 数据 集 具 有 
较 高 的 准确 性 。Kim 等 "在 CMORPH 降雨 数据 集 
的 基础 上 ,采用 El 模型 对 美国 地 区 的 降雨 侵蚀 力 
进行 估算 ,得 出 CMORPH 对 沿海 地 区 的 降雨 侵蚀 力 
存在 低估 现象 ,对 于 站 点 稀 玻 的 地 区 能 够 较 好 地 反 
映 其 降雨 侵蚀 状况 。 

综 上 所 述 , 已 有 的 降雨 侵蚀 力 研 究 主 要 存在 以 
下 2 个 问题 :(1) 使 用 气象 站 数据 计算 降雨 侵蚀 力 
仅 能 反映 站 点 附近 有 限 区 域 的 侵蚀 状况 ,而 在 进行 
大 尺度 计算 时 存在 空间 分 异性 的 问题 ”。 卫 星 降 
雨 资料 具有 空间 覆盖 范围 广 ,捕捉 降水 连续 时 空 变 
化 的 特性 ,特别 是 在 站 点 稀 玻 区 具有 明显 的 优势 ， 
但 它 作 为 一 种 间接 测 雨 方法 , 受 传 感 澡 性能. 反 演 
算法 等 限制 ,在 测量 过 程 中 易 产 生 误差 ,估算 的 降 
十 侵蚀 力 精度 有 限 “”,(2) 受降 雨 资料 获取 的 限制 ， 
前 人 对 于 降雨 侵 刨 力 的 研究 多 是 在 局 地 尺度 进行 ， 
且 缺 乏 对 中 国 近 20 a 来 长 序列 降 南 侵蚀 力 的 深入 分 
析 。 在 少 有 的 针对 全 国 尺度 的 研究 中 ,多 是 采用 站 
点 估算 和 搬 值 结合 获得 降雨 侵蚀 力 的 空间 分 布 ,未 
考虑 站 点 稀 玻 区 降雨 侵蚀 力 估算 的 问题 ,而 在 使 用 
卫星 数据 的 相关 研究 中 则 无 法 避免 该 资料 精度 低 
的 缺陷 。 

基于 以 上 存在 的 问题 ,本 次 研究 使 用 国家 气象 
言 息 中 心 研 制 的 CLDAS 多 源 融合 降水 数据 集 , 通 过 
Pa 模型 对 中 国 的 降雨 侵蚀 力 进行 估算 ,从 不 同 的 
时 间 空间 尺度 评估 CLDAS 降雨 侵蚀 力 的 整体 精 
BE ,并 结合 降雨 量 侵蚀 性 降雨 次 数 和 侵蚀 密度 ,对 
研究 区 的 降雨 侵蚀 力 特征 进行 分 析 ,为 中 国 区 域 的 
水 土 保 持 研 究 以 及 实践 工作 提供 科学 参考 依据 。 


1 数据 与 方法 


1.1 数据 来 源 

1.1.1 CLDAS 多 源 融合 降水 ”国家 气象 信息 中 心 
Sun 等 中 使 用 多 重 网 格 变 分 分 析 方 法 中 ,以 
CMORPH 卫星 降水 资料 和 美国 航空 航天 局 研发 的 


MERR2(Modern-Era retrospective analysis for research 
and applications, Version 2) 再 分 析 降 水 资料 为 联合 
背景 场 ,融合 2400 多 个 国家 级 站 点 观测 数据 和 6x 
10' 多 个 区 域 自动 站 的 降水 数据 制作 了 一 套 空间 分 
PERI 6.25 kmx6.25 km, 时 间 分 辩 率 为 1 的 多 源 
融合 降水 产品 CLDAS-prep (CLDAS Multi-source pre- 
cipitation fusion dataset over China ) ,该 数据 可 通过 国 
家 气象 科学 数据 中 心 (http://data.cma.cn/) 申 请 下 载 ， 
本 研究 所 用 资料 的 时 序 长 度 为 2001 一 2020 年 。 
11.2 地 面 气 象 站 数据 地 面 气象 站 资料 来 自 国家 
气象 信息 中 心 2001 一 2020 年 中 国 国家 级 地 面 气象 
站 逐 小 时 降水 数据 ,可 从 国家 气象 科学 数据 中 心 
(http://data.cma.cn/) 申 请 下 载 。 
1.2 研究 方法 
1.2.1 降雨 侵蚀 力 计算 ”降雨 侵蚀 力 计算 采用 Wis- 
chmeier 经 典 算 法 公式 ": 
R= m$ Zen) a) 

式 中 :及 为 年 平均 降雨 侵蚀 力 因子 (MJ'mm'hm 2?.h…， 
a );n 为 计算 R 值 的 总 年 数 ;mi 为 第 i 年 内 的 侵蚀 性 
降雨 次 数 ;为 次 降雨 总 动能 (MJ:hnr) ;J 为 次 降雨 
最 大 30 min 雨 强 (mm*h'); 文 草 中 的 /表示 最 大 nn 
分 钟 雨 强 (mm'h- )。 

由 于 CLDAS 融 合 降水 的 时 间 分 辩 率 为 60 min, 
故 本 研究 参考 Yue 等 "对 小 时 降水 提出 的 改进 公式 : 


R= mÈ Zer) xl1871 (2) 


i=1| j=l 


E=% ep, (3) 
k=l 


SUP : m AGU HE BEB PE m BS Be faz AE BE 
雨 的 划分 标准 是 单 次 降雨 过 程 的 降雨 量 达 到 12.70 
mm) ;e 为 第 上 个 阶段 的 单位 降雨 动能 (CMJ hm 
mm) ;pi 为 第 个 阶段 的 降雨 量 (mm )。 

每 个 阶段 的 单位 降雨 动能 (e) 的 计算 公式 如 下 

e =0.29(1 - 0.72e °°”) (4) 

式 中 :i 为 每 个 阶段 内 的 降雨 强度 (mm:h)。 
1.2.2 侵蚀 密度 ”侵蚀 密度 (ED,MJ.hm”…h') 是 降 
雨 侵蚀 力 与 降雨 量 的 比值 ,反映 的 是 单位 雨量 的 侵 
刨 强度。 


ED=R/P (5) 
式 中 :R 为 降雨 侵 刨 力 值 (MJ.mm:hm*…h');P 为 降 


雨量 (mm )。 
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1.2.3 降雨 侵蚀 力 评 估 根据 气象 观测 站 的 经 纬度 
(图 1) ,将 各 站 点 匹配 至 CLDAS 降雨 侵蚀 力 格 点 上 ， 
采用 双 线 性 插值 提取 出 降雨 侵蚀 力 值 ,利用 皮尔 进 
相关 系数 (7) .相对 误差 (BIAS) 指 标 评估 CLDAS 降 
十 侵蚀 力 数据 的 精度 。 利 用 皮尔 人 逊 相关 系数 捅 述 
降雨 侵蚀 力 值 与 真 值 之 间 的 吻合 程度 ;利用 相对 误 
差 描述 降雨 侵蚀 力 的 误差 特征 。 具 体 公 式 如 下 : 


r= et (6) 


x 100% (7) 


式 中 :NN 为 样本 数目 ; C; 为 CLDAS 降雨 侵蚀 力 值 ; C 
为 CLDAS 降雨 侵蚀 力 的 平均 值 ; G, 为 实测 降雨 侵 
蚀 力 值 ; G 为 站 点 的 实测 降雨 侵蚀 力 平均 值 。 


2 结果 与 分 析 


2.1 降雨 侵蚀 力 计 算 结 果 分 析 
2.1.1 降雨 侵蚀 力 精度 分 析 图 2a~d 为 CLDAS 降 
雨 侵蚀 力 和 站 点 侵蚀 力 在 季节 尺度 的 散 点 图 ,根据 


。 气象 站 点 
国界 
- - 未 定 国界 
一 省 界 
台湾 省 资料 暂 缺 
0 500 km 
= 人 sh 


HAS: GS(2022)3443 号 


二 者 的 一 元 线性 回归 ,CLDAS 降雨 侵蚀 力 整体 上 略 
低 于 地 面 站 数据 ,但 其 相关 系数 均 达 到 0.80 以上， 
表现 出 良好 的 线性 关系 , 春 、 夏 、 秋 季 的 相关 系数 达 
到 0.90 以 上 ,而 冬季 由 于 固态 降水 观测 的 误差 大 ， 
其 相关 系数 低 于 其 他 季节 ,为 0.82。 从 相对 误差 来 
看 ,春季 的 误差 最 低 ,为 24.67% ,冬季 的 误差 最 大 ， 
为 33.40% ,其 次 是 夏季 和 秋季 ,分 别 为 29.71% 和 
32.30%。 总 体 上 ,在 降雨 丰富 的 季节 ,CLDAS 降 雨 
侵蚀 力 的 适应 性 较 好 ,可 以 准确 呈现 出 降雨 侵蚀 力 
的 季节 性 变化 。 

图 2e~f 为 CLDAS 降雨 侵蚀 力 和 站 点 侵蚀 力 在 
月 尺度 .年 尺度 的 一 元 线性 回归 ,其 中 月 尺度 的 相 
关系 数 达 到 了 0.92 ,相对 误差 为 29.85% ;年 尺度 的 
相关 系数 为 0.96, 相 对 误差 为 28.80%。 由 于 降水 数 
据 在 月 尺度 上 的 随机 误差 普遍 高 于 年 尺度 ,同时 误 
差 会 随 着 模型 传递 ,从 而 影响 最 终 的 降雨 侵蚀 力 结 
He ,因此 年 尺度 的 降 十 侵蚀 力 精度 总 体 优 于 月 
尺度 。 

Kim 等 的 利用 CMORPH 卫星 降水 产品 对 美国 
地 区 的 降雨 侵蚀 力 进 行 评估 ,结果 表明 CMORPH 4E 
降雨 侵蚀 力 的 相关 系数 为 0.80, 相 对 误差 为 
64.00% , CLDAS 降水 在 CMORPH 资料 的 基础 上 经 
过 站 点 融合 后 ,使 得 降雨 侵蚀 力 的 相对 误差 显著 降 


图 1 中 国 气象 站 点 分 布 


Fig. 1 Distribution of meteorological stations in China 
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图 2 站 点 降雨 侵蚀 力 与 CLDAS 降雨 侵蚀 力 在 春 . 夏 、 秋 、 冬 季 以 及 月 尺度 \ 年 尺度 的 线性 回归 


Fig.2 Linear regressions of weather station rainfall erosivity and CLDAS rainfall erosivity 


in spring, summer, autumn and winter, as well as monthly and annual scales, respectively 


低 , 相 关系 数 从 0.80 提 升 至 0.96。 

综 上 ,CLDAS 降雨 侵蚀 力 在 不 同 的 时 间 尺 度 上 

均 有 较 好 的 精度 ,能 够 克服 气象 站 和 卫星 降雨 资料 
各 自 存在 的 不 足 , 对 于 揭示 降雨 侵蚀 力 的 时 空 变化 
和 精确 佑 算 土 壤 侵蚀 等 方面 具有 显著 优势 。 
2.1.2 年 降雨 侵蚀 力 分 析 图 3a~c 为 年 降雨 侵蚀 
J 年 降雨 量 年 侵蚀 性 降雨 次 数 在 2001—2020 年 
的 变化 趋势 。 由 于 中 国 不 同 地 区 的 地 形 和 和 气候 存 
在 差异 ,使 得 各 区 域 的 降雨 量 和 降雨 强度 等 表现 不 
同 ,本文 以 年 降雨 量 <400.00 mm, 400.00~800.00 
mm ,>800.00 mm 为 标准 ,将 我 国 划 分 为 低 、 中 、 高 3 
个 不 同 的 雨量 分 区 ,分 析 在 不 同 的 雨量 条 件 下 各 指 
标的 变化 情况 。 


总 体 上 3 个 雨量 区 的 年 降雨 侵蚀 力 年 降雨 量 
和 年 侵蚀 性 降雨 次 数 的 变化 趋势 基本 一 致 。 在 低 
雨量 区 (<400.00 mm) ,年 降雨 侵蚀 力 为 98~549 MJ- 
mm:hm?:h!+a!, 4E KER a AY 162.00~225.00 mm ,年 
侵蚀 性 降雨 次 数 为 3~6 次 ;在 中 雨量 区 (400.00~ 
800.00 mm) ,年 降雨 侵蚀 力 为 649~947 MJ - mm: 
hm”h"'…a', 年 降雨 量 为 536.72~579.84 mm ,年 侵蚀 
性 降雨 次 数 为 10~15 次 。 由 图 3a~c 可 知 ,3 个 指标 
在 低 雨 量 区 和 中 雨量 区 的 变化 幅度 基本 稳定 ,而 在 
高 雨量 区 的 年 际 波动 剧烈 。 

在 高 雨量 区 ,降雨 侵蚀 力 最 低 为 2853 MJ-mm- 
hm?-+h7!+a7!(2003 Œ) ,最 高 达到 4250 MJ-mm-hm?: 
ha (2008 年 )。 同 时 ,2003 年 的 降雨 量 和 侵蚀 性 
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图 3 2001 一 2020 年 的 降雨 侵蚀 力 、 降 雨量 侵蚀 性 降雨 次 数 和 暴雨 量 在 不 同 雨 量 区 的 逐年 变化 


Fig. 3 Annual variation of rainfall erosivity, rainfall, erosive rainfall events and heavy 


rainfall in different rainfall zones from 2001 to 2020 


降雨 次 数 最 低 ,分别 为 1211.54 mm Fil 26 YK ; 2002 4F 
的 降雨 量 最 高 ,达到 1529.67 mm; 2020 年 的 侵蚀 性 
降雨 次 数 最 多 ,为 41 次 。 

高 雨量 区 的 部 分 年 份 呈 现 出 低 雨 量 .高 侵蚀 力 
的 特征 ,如 2020 年 和 2016 年 的 年 降雨 侵蚀 力 分 别 
为 4225 Mj:mm:hm ha!.4094 MJ.mm.hm2.h 
a', 但 2020 年 的 年 均 降雨 量 为 1329.15 mm, 比 2016 
年 减少 了 180.00 mm。 根 据 图 3c 可 得 ,2020 年 的 侵 
蚀 性 降雨 发 生 最 为 频繁 ,结合 侵蚀 力 计算 公式 ,其 


严重 影响 。 然 而 从 图 3a 和 图 3c 可 知 ,2019 年 的 侵 
蚀 性 降雨 事件 较 少 ,同时 该 年 的 降雨 侵蚀 力 整体 低 
于 其 他 年 份 ,未 出 现 类 似 2020 年 的 高 暴雨 伴随 高 侵 
蚀 力 的 现象 。 由 此 可 见 , 降 雨量 一 定时 ,侵蚀 性 降 
雨 次 数 和 暴雨 过 程 协同 决定 降雨 侵蚀 力 的 大 小 。 
2001 一 2020 年 的 多 年 平均 降雨 侵蚀 力 分 布 特 
征 如 图 4 所 示 , 降 十 侵蚀 力 分 布 图 的 空间 分 辩 率 为 
6.25 kmx6.25 km ,平均 降雨 侵蚀 力 为 1238 MJ-mm- 
hm?+h! -a~ ,标准 差 为 1802 MJ.mm.hn .ha , KE 


中 60 min 最 大 雨 强 对 降雨 侵蚀 力 的 影响 权重 最 高 ， 
说 明 2020 年 的 降雨 特点 为 历时 短 .强度 大 。 

图 3d 为 2001 一 2020 年 的 逐年 暴雨 总 量 (24n 累 
积 降水 超过 50.00 mm 的 降雨 事件 ) ,其 中 2020 年 的 
暴雨 量 明 显 高 于 2016 年 ,可 印证 前 文 2020 年 强 降 
雨 事件 密集 的 结论 。 从 全 国 尺度 来 看 ,2020 年 的 暴 
雨量 最 多 ,而 在 高 雨量 区 ,2019 年 的 暴雨 趋势 呈 单 
峰 状 ,主要 原因 是 2019 年 受 人 台风 “ 利 奇 马 "影响 ,长 
江 .黄河 等 流域 迎 来 最 强暴 十 ,降雨 总 量 达 1336.9x 
10m ,对 山东 省 2 浙江 省 等 地 区 的 水 土 流失 造成 


雨 侵蚀 力 随 空间 位 置 而 变化 。 总 体 上 看 ,中国 西 部 
的 降雨 侵蚀 力 较 低 , 东 南 沿海 地 区 的 降雨 侵蚀 力 较 
高 。 内 蒙古 自治 区 .青海 省 .新 疆 维吾尔 自治 区 西 
藏 自治 区 等 内 陆地 区 的 降雨 侵蚀 力 低 于 500 MJ - 
mm*hm”…h'…a', 而 降雨 侵蚀 力 高 于 5000 MJ mm- 
hm? +h! +a ! 的 地 区 主要 分 布 在 广西 广东 ,福建 等 沿 
海 省 份 。 我 国 的 东北 地 区 .中 部 地 区 和 西南 地 区 的 
降雨 侵蚀 力 大 部 分 为 1000~5000 MJ *mm*hm™®- h“. 
ao AR 1 可 知 我 国 各 省 份 ( 地 区 ) 的 年 均 侵 蚀 力 及 
其 变化 范围 ,台湾 省 .香港 .澳门 海南 省 显著 大 于 
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图 4 2001 一 2020 年 中 国 年 降雨 侵蚀 力 空间 分 布 
Fig.4 Spatial distribution of annual rainfall erosivity in China from 2001 to 2020 


#1 2001 一 2020 年 中 国 各 省 份 ( 地 区 ) 年 降雨 侵蚀 力 
Tab. 1 Annual rainfall erosivity in each province (region) of China from 2001 to 2020 


省 份 (地 区 ) E E r PEN ie A es. 
范围 平均 值 范围 平均 值 
北京 市 334~2777 1552 湖南 省 1529~6013 3405 
天 津 市 1517~3232 2042 广东 省 2573~15716 7009 
河北 省 140~3100 1153 广西 省 2003~16295 5331 
山西 省 211~1805 620 海南 省 2947~12717 8710 
内 蒙古 1~1642 348 重庆 市 1424~4811 2616 
辽宁 省 383~4121 1615 四 川 省 73~5678 1488 
吉林 省 391~2071 859 贵州 省 678~5935 27908 
黑龙 江 省 223~2198 834 云南 省 159~6715 1674 
上 海 市 2423~4399 3371 西藏 2~53718 583 
江苏 省 2085~5533 3490 陕西 省 229~3599 878 
浙江 省 1642~7449 4129 甘肃 省 0~2314 228 
安徽 省 1802~8002 3549 青海 省 0~16551 309 
福建 省 2597~8529 4713 宁夏 60~1016 251 
江西 省 2728~7664 4587 新 疆 0~10163 63 
山东 省 1307~4077 2382 台湾 3207~14078 7966 
河南 省 494~4030 1840 香港 2928~9190 5771 
湖北 省 719~5702 2763 澳门 6364~12119 8747 


其 他 省 市 ,平均 侵蚀 力 为 5000~12000 MJ:mm: hm mm 等 降水 量 线 划 分 ,此 线 以 西 以 北 的 地 区 深 处 内 
ha ,其 次 是 广东 省 ,福建 省 等 沿海 省 份 。 陆 , 降 水 稀少 ,因此 侵 刨 力 较 小 , 均 小 于 500 MJ- 
综 上 ,年 降雨 侵蚀 力 的 空间 分 布 大 致 以 400.00 mm.hm2.h.a ,同时 年 际 间 的 变化 相对 稳定 ;此 线 
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以 东 以 南 的 地 区 ,年 降雨 侵蚀 力 较 大 , 且 变 率 也 大 。 

2.1.3 月 降雨 侵蚀 力 分 析 图 5、 图 6 分 别 为 2001 一 
2020 年 1 一 6 月 和 7 一 12 月 的 逐 月 降雨 侵蚀 力 。 冬 
季 (12.1. 2 月) , 除 沿海 省 份 外 ,中 国 的 大 部 分 地 区 
以 固态 降水 为 主 ,发 生 强 降雨 的 频次 很 低 ,所 以 冬季 
大 范围 内 的 降雨 侵蚀 力 小 于 25 MJ -mmhmmh .an'; 
但 在 气候 较为 湿润 的 东南 沿海 地 区 (浙江 省 .福建 
省 等 ) ,侵蚀 性 降雨 的 出 现 导致 相对 较 高 的 降雨 侵 
刨 力 。 春 季 (3、4、5 月 ) ,我 国 华南 地 区 在 副热带 高 


By 


压 疹 雨 带 影 响 下 进入 雨季 ,发 生 强 降雨 的 频次 增 
多 ,降雨 侵蚀 力也 逐渐 增 大 ,如 广东 省 .广西 省 等 地 
区 ,降雨 侵蚀 力 相 较 于 冬季 有 明显 的 增强 趋势 。 

随 着 副 高 压 辐 散 的 暧 湿 气 流 不 断 问 北 移动 ,我 
国 长 江 中 下 游 地 区 5 月 底 一 6 月 初 进入 雨季 ,降雨 的 
强度 和 范围 也 大 幅 增加 ,如 湖南 省 .湖北 省 .安徽 省 
等 地 区 ,降雨 侵蚀 力 大 于 650 MJ-mm-hm?-h'+a', 
北 太 平 洋 副热带 高 压 的 不 断 北 移 导 致 锋面 区 域 形 
成 的 降雨 带 向 北 移动 ,直到 7.8 月 在 黄河 .淮河 流域 


加 
-未定 国界 SO Ss N 


图 5 2001 一 2020 年 1 一 6 月 中 国 降雨 侵蚀 力 的 空间 分 布 


Fig. 5 Spatial distributions of rainfall erosivity in China in January-June from 2001 to 2020 


202210.00194v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


(e) 11 月 


(d) 10 月 


图 例 - 月 降雨 侵 亿 力 
> OS 8S 8 G8 OP 8 /M3-mm:hm?:hr!-a! 


SMS KS 7 


图 6 2001 一 2020 年 7 一 12 月 中 国 降雨 侵蚀 力 的 空间 分 布 
Fig.6 Spatial distributions of rainfall erosivity in China in January-June from 2001 to 2020 


停滞 ,造成 多 次 的 ,大 范围 的 强 降雨 。 进 入 华北 雨 
季 , 月 降雨 侵蚀 力 急剧 增加 ,尤其 是 在 东北 、 华 北 以 
及 华中 的 部 分 地 区 ,如 山西 省 .山东 省 ` 河 北 省 等 地 
区 的 降雨 侵蚀 力 超过 450 MJ:mm:hm?:h'.a'。10 
月 , 随 着 副热带 高 压 的 迅速 南 退 ,我国 雨季 随 之 结 
束 , 发 生 强 降雨 事件 的 频次 大 大 减少 ,此 时 我 国 以 
西 以 北 降雨 侵 刨 力 大 部 分 小 于 25 MJ-mm-hm®?-h'- 
a', 以 东 以 南 的 地 区 降雨 侵蚀 力 小 于 150 MJ mm 
hm ha, 总体 上 降雨 侵蚀 力 的 空间 分 布 仍 呈 现 
东南 高 .西北 低 的 特点 。 


2.2 侵蚀 密度 结果 分 析 
2.2.1 年 侵蚀 密度 分 析 2001 一 2020 年 的 侵蚀 密度 
如 图 7a 所 示 。 中 国 的 年 侵蚀 密度 值 为 1.40 MJ-hm?- 
pa, 东南 沿海 地 区 部 分 省 份 的 年 侵蚀 密度 达到 
3.50 MJ .hn .ha 以 上 ,如 广东 省 海南 省 .台湾 
省 等 地 区 。 华 中 和 华东 地 区 非 沿 海 省 份 的 侵蚀 密 
度 多 稳定 在 2.30~5.00 MJ -hm?-h'+ a! Z fia] , RAL 
区 的 侵蚀 密度 为 0.45~1.50 MJ.hm2.h -ao 

表 2 为 各 省 份 ( 地 区 ) 的 年 均 降 雨量 、 年 均 侵蚀 
密度 和 年 均 暴 雨量 。 海 南 省 的 年 均 侵 蚀 密度 .年 均 
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图 7 2001 一 2020 年 中 国 年 侵蚀 密度 空间 分 布 
Fig.7 Spatial distribution of annual erosivity density in China from 2001 to 2020 


#2 2001 一 2020 年 中 国 各 省 份 (地 区 ) 年 均 降雨 量 ,年 均 侵 蚀 密度 ,年 均 暴雨 量 
Tab. 2 Average annual rainfall, average annual erosion density, and average annual 
heavy rainfall by province (region) in China from 2001 to 2020 


北京 市 553.58 2.75 75.19 湖南 省 1452.26 2.33 233.48 
天 津 市 556.38 3.67 93.98 广东 省 1889.15 3.66 503.05 
河北 省 496.50 2.20 52.62 广西 省 1613.90 3.26 361.97 
山西 省 498.18 1.22 29.05 海南 省 1910.84 4.55 619.89 
内 蒙古 294.58 1.06 11.88 重庆 市 1197.47 2.18 179.20 
辽宁 省 653.05 2.41 107.45 四 川 省 900.14 1.47 88.83 
吉林 省 600.07 1.46 43.84 贵州 省 1169.93 2.33 158.69 
黑龙 江 省 550.93 1.50 34.57 云南 省 1092.31 1.47 81.70 
上 海 市 1305.19 2.58 231.82 西藏 388.74 1.08 29.14 
江苏 省 1112.86 3.18 237.27 陕西 省 651.67 1.32 53.28 
浙江 省 1686.28 2.44 297.65 WE 287.51 0.63 7.61 
安徽 省 1266.17 2.82 261.37 青海 省 329.23 0.84 8.62 
福建 省 1764.18 2.68 352.82 宁夏 291.60 0.86 8.48 
江西 省 1734.80 2.63 327.61 新 疆 140.19 0.47 1.83 
山东 省 704.64 3.37 147.78 台湾 1709.28 4.62 711.58 
河南 省 764.97 2.34 123.28 香港 1624.13 3.51 468.18 
湖北 省 1155.33 2.34 191.13 澳门 1943.85 4.36 538.17 
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降雨 量 和 年 均 暴雨 量 最 高 ,分 别 为 4.55 MJ .hm ”.h 
a',1910.84 mm 和 619.89 mm ,新 疆 维 吾 尔 自治 区 的 
年 均 侵蚀 密度 最 低 , 为 0.47 MJ hm? +h! +a’, [ej EE 
年 均 降 雨量 也 最 低 。 贵 州 省 和 云南 省 的 年 均 降雨 
量 差异 小 ,分 别 为 1169.93 mm 和 1092.31 mm ,年 均 
侵蚀 密度 为 2.33 MJ- hm? +h! +a Fl 1.47 MJ-hm?-h 
a ,前 者 的 降雨 量 是 后 者 的 1.07 倍 ,但 年 均 侵 蚀 密 
度 是 1.58 倍 ,由 此 可 得 出 ,贵州 省 的 单位 降雨 量 能 
带 来 较 高 的 降雨 侵蚀 力 。 从 两 省 的 暴雨 量 可 知 , 贵 
州 省 的 年 均 暴 雨量 为 158.69 mm ,云南 省 的 年 均 暴 
雨量 为 81.7 mm ,两 者 的 比例 为 1.93:1.00 ,说 明年 均 
暴雨 量 对 年 均 侵 蚀 密 度 的 影响 更 敏感 。 结 合 侵蚀 
密度 图 (图 7) 分 析 , 贵 州 省 的 侵蚀 密度 更 高 .范围 更 
广 , 因 此 降雨 侵蚀 力 更 大 ,具有 这 种 关系 的 省 份 还 
有 辽宁 省 和 吉林 省 ,两 省 的 年 均 降 雨量 分 别 为 
653.05 mm 和 600.07 mm, 但 年 均 侵 蚀 密 度 差异 显 
著 , 分 别 为 2.41 MJ.hm2.pa 和 1.46 MJ.hm2.h-， 
a ,而 对 应 的 年 均 暴雨 量 分 别 为 107.45 mm 和 43.84 
mm。 综 上 可 知 , 暴 雨量 越 高 ,侵蚀 密度 越 大 。 


(a) 春季 


图 例 

一 国界 四 

--- 未 定 国界 © sss 
审 图 号 : GS(2022)3443 号 一 省 界 Sr > 


(b) 夏季 


2.2.2 季节 侵蚀 密度 分 析 ”由 季节 侵蚀 密度 图 (图 
8) 可 知 ,不同 季 节 的 侵蚀 密度 空间 分 布 差异 较 大 。 
春季 ,东南 沿海 省 份 以 及 湖南 省 贵州 省 .江西 省 的 
侵蚀 密度 大 于 1.80 MJ-hm? hla! ,江苏 省 安徽 省 
等 地 区 的 侵蚀 密度 在 1.20~1.80 MJ -hmh -a 2M] 
浮动 ,全 国 的 平均 侵蚀 密度 为 0.85 MJ .hm h-a, 
标准 差 为 1.44 MJ-hm?-h'+a!, B22, (2 phe 
值 区 迅速 扩张 ,平均 侵蚀 密度 为 1.81 MJ.hm ha ， 
标准 差 为 1.19 MJ.hna2.ha ,华东 、 华 中 以 及 东北 
地 区 的 侵蚀 密度 高 于 2.30 MJ.hm: .ha'。 秋 季 ， 
四 川 省 东部 和 东南 沿海 地 区 的 侵蚀 密度 为 2.30~5 
MJ .hn2.ha ,海南 省 和 台湾 省 的 局 部 地 区 超过 
5.00 MJ"hm .ha ,全 国平 均 侵蚀 密度 为 0.88 MJ- 
hm-2 .ha ,标准 差 为 0.92 MJ-hm?-ht-a!, 冬季， 
由 于 西伯 利 亚 冷 高 压 控制 ,全国 的 侵蚀 密度 整体 下 
降 , 低 于 1.50 MJ.hm2.h 2a ,平均 侵蚀 密度 为 1.10 
MJ .hoa-2.ha ,标准 差 为 1.27 MJ-hm?-h'-a', Re- 
nard 等 "研究 发 现 , 季 节 侵 刨 密度 与 30 min 两 强 成 
正比 ,可 反应 降雨 强度 的 季节 性 变化 ,总 体 上 夏季 


es 


: -2.h-l.a-! 
$ /MJ'hm ha 


图 8 2001 一 2020 年 中 国 季节 侵蚀 密度 空间 分 布 
Fig. 8 Spatial distributions of seasonal erosivity density in China from 2001 to 2020 
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的 侵蚀 密度 最 高 ,春季 和 秋季 次 之 ,冬季 最 低 。 

当 侵蚀 密度 大 于 1 MJ hm? +h! +a Ht, PAE ERD 
量 会 带 来 较 高 的 降雨 侵蚀 力 ,高 侵蚀 性 的 降水 冲击 
长 期 干燥 的 土壤 ,通常 会 造成 巨大 的 损害 ,同时 有 
研究 表明 ,降水 集中 程度 增 大 会 增加 干旱 发 生 的 可 
能 性 2 ,因此 本 研究 对 于 降雨 量 稀 少 但 侵蚀 密度 高 
且 植 被 覆盖 度 低地 区 的 土壤 侵蚀 防护 工作 具有 指 
示意 义 。 


3 讨论 


气象 站 观测 目前 仍 是 估算 降雨 侵蚀 力 的 主流 
数据 ,然而 站 点 资料 虽 精 度 高 但 不 适用 于 大 区 域 范 
围 内 的 降雨 侵蚀 力 估 算 , 它 仅 能 反映 站 点 附近 区 域 
的 降水 现状 ,卫星 降雨 产品 虽 能 有 效 再 现 降水 的 时 
空 分 布 , 但 受 反 演算 法 和 传感器 的 影响 ,其 精度 一 
直 是 在 各 专业 领域 应 用 时 的 瓶 贷 。 近 10 a 来 ,多 源 
降水 资料 融合 的 理念 为 研发 高 精度 的 降水 产品 提 
供 了 新 的 思路 ,融合 降水 数据 是 通过 一 定 的 融合 方 
法 ,在 同一 时 空 尺度 上 将 站 点 观测 数据 、 卫 星 降水 
产品 、 雷 达 降 水 产品 等 不 同 来 源 的 降水 资料 进行 融 
合 , 融 合 后 的 降水 具有 多 种 资料 的 优势 ,更 加 接近 
降水 的 真实 空间 分 布 沾 。CLDAS 多 源 融 合 降水 数 
据 由 国家 气象 信息 中 心 自主 研制 ,被 国内 外 学 者 广 
泛 应 用 在 陆 面 模式 习 、 山 洪 凶 以 及 干旱 监测 轧 等 领 
域 。 本 研究 首次 将 CLDAS 融合 降水 应 用 在 降雨 侵 
刨 力 的 研究 中 ,并 首次 利用 逐 小 时 国家 站 数据 评估 
CLDAS 降雨 侵蚀 力 的 精度 ,为 该 数据 集 在 土壤 侵蚀 
中 的 应 用 提供 参考 依据 。 

估算 降雨 侵蚀 力 的 模型 主要 分 为 雨 强 模型 和 
简易 模型 '*, 雨 强 模型 以 Wischmeier 经典 算法 模 
型 “为 代表 ,Wischmeier 等 “利用 美国 的 35 个 水 土 
保持 站 、8250 个 小 区 的 降雨 和 侵蚀 实测 资料 ,以 次 
降雨 为 单位 ,将 降雨 总 动能 (5) 与 30 min 最 大 雨 强 
(Too) HY HE AR Bl 作为 降雨 侵 刨 力 指 标 ,此 方法 一 直 
被 各 国学 者 作为 标准 法 沿用 ,然而 它 一 般 需 要 高 时 
间 分 状 率 的 降雨 过 程 资 料 , 且 费时 费力 ,对 于 长 时 
间 序 列 或 大 尺度 区 域 范 围 内 的 降雨 侵蚀 力 计算 难 
以 实现 。 由 于 日 尺度 、 月 尺度 降雨 资料 的 易 获取 
性 , 章 文 波 等 "\ 王 万 中 等 "建立 了 简易 估算 模型 ， 
例如 日 雨量 模型 和 月 雨量 模型 ,但 这 些 模 型 大 多 是 
在 雨 强 模型 的 基础 上 通过 统计 拟 合 出 最 优 参数 , 往 


往 在 研究 区 内 效果 良好 ,而 在 研究 区 之 外 的 范围 适 
用 性 较 低 ,例如 日 雨量 模型 在 降雨 量 丰 富 的 地 区 ， 
计算 得 到 的 降雨 侵蚀 力 误差 相对 小 ;在 降雨 量 稀少 
的 地 区 ,模型 计算 的 侵蚀 力 误差 较 大 史 ,因此 将 简 
易 模 型 应 用 在 全 国 尺度 存 在 一 定 的 局 限 性 裤 。 本 
人 研究 依靠 国家 气象 信息 中 心 高 性 能 服务 器 ,基于 py- 
thon 开发 了 适用 于 CLDAS 融合 降水 产品 的 Kw 模 
型 ,克服 了 雨 强 模型 获取 数据 难 和 计算 费时 的 问 
题 ,并 能 够 在 全 国 尺度 展开 长 序列 研究 。 

由 于 降雨 侵蚀 力 逐 年 的 时 空 变化 显著 ,至少 需 
要 15 a 以 上 的 降水 资料 才能 获取 具有 代表 性 的 年 降 
雨 侵 刨 力 结果 ”。 受 降雨 资料 和 计算 资源 的 限制 ， 
以 往 对 于 全 国 尺度 的 降雨 侵蚀 研究 仅 分 析 单个 年 
份 的 状况 ,难以 在 全 国 尺 度 展开 长 序列 的 降雨 侵蚀 
力 分 析 , 本 文 使 用 20 a 的 融合 降水 资料 估算 中 国 的 
降雨 侵蚀 力 , 在 长 时 间 尺 度 上 的 结果 更 具有 代 
表 性 。 

此 外 ,本 文 的 研究 还 存在 一 些 不 足 , 多 数 卫 星 
降水 产品 可 能 高 估 湖 泊 和 水 库 区 域 的 降水 ,由 于 
CLDAS 融 合 降水 以 CMORPH 为 背景 场 ,因此 西部 湖 
泊 区 域 的 降雨 侵蚀 力 存在 明显 的 高 估 现 象 , 比较 典 
型 的 有 纳木错 和 色 林 错 湖泊 等 。 张 黎明 ”等 利用 
1997 一 2003 年 的 自 建 气象 站 降雨 资料 ,研究 了 不 同 
降雨 强度 对 估算 侵蚀 力 的 影响 ,结果 表明 7o>7Fs>Po> 
J。 因此 CLDAS 融合 降水 的 时 间 分 辩 率 仍 有 所 不 
足 , 未 来 可 利用 分 钟 级 的 气象 站 观测 资料 对 CLDAS 
融合 降水 进行 时 间 降 尺度 ,提升 其 在 降雨 侵蚀 力 
估算 中 的 精度 。 


4 结 论 

针对 气象 站 和 卫星 降雨 资料 在 降雨 侵蚀 力 应 
用 中 存在 无 法 反映 空间 异 质 性 和 数据 精度 差 的 问 
题 ,本 文 使 用 CLDAS 多 源 融 合 降水 数据 结合 地 面 气 
象 站 数据 ,采用 Els 模型 ,对 中 国 的 降雨 侵蚀 力 进 行 
分 析 评 佑 ,得 出 以 下 结论 : 

(1) CLDAS 融 合 降水 估算 的 降雨 侵蚀 力 在 全 国 
范围 内 具有 较 好 的 适应 性 ,与 地 面 实测 结果 在 不 同 
的 时 间 尺 度 具 有 和 良好 的 线性 回归 关系 ,相关 系数 均 
达到 0.80 以 上 ,相对 误差 显著 降低 ,可 应 用 于 长 序 
列 大 尺度 范围 内 的 降雨 侵蚀 力 研究 。 

(2) 2001 一 2020 年 ,不同 雨量 区 的 年 降雨 侵蚀 
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力 年 降雨 量 \ 年 侵蚀 性 降雨 次 数 的 变化 趋势 基本 
一 致 。 在 低 雨 量 区 和 中 雨量 区 ,3 个 指标 逐年 变化 
稳定 ,而 在 高 雨量 区 的 年 际 波动 剧烈 ,侵蚀 性 降雨 
次 数 和 暴雨 过 程 协 同 影响 降雨 侵蚀 力 的 大 小 。 

(3) 中 国 的 降雨 侵蚀 力 在 空间 上 大 致 呈现 出 东 
南 沿海 地 区 高 西北 内 陆地 区 低 的 特征 。 时 空 上 ， 
全 国 的 侵蚀 性 降雨 集中 在 5 一 8 月 ,各 月 之 间 的 侵蚀 
力 变 化 幅度 大 ,年 均 降 十 侵蚀 力 为 1238 MJ- mm: 
hm? +h'+a'; L4 400.00 mm 等 降水 量 线 划 分 ,此 线 以 
西 以 北 降水 稀少 ,侵蚀 力 低 , 均 小 于 500 MJ -mm 
hm- h” -a ,此 线 以 东 以 南 的 地 区 受 季 风 影 响 ,降雨 
量 增多 ,侵蚀 力 较 高 。 

(4) 通过 对 年 降雨 量 .年 侵蚀 密度 和 年 暴雨 量 
进行 分 区 定量 分 析 , 得 出 暴雨 量 与 侵蚀 密度 成 正 相 
关 即 年 降雨 量 一 致 时 ,暴雨 量 越 高 ,侵蚀 密度 越 大 。 
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Abstract: Rain gauge data is still the most popular data for estimating rainfall erosivity. However, although this 
precipitation data is highly accurate, it is not suitable for estimating rainfall erosivity on a large, regional scale. It 
can only reflect the situation of rainfall erosivity in the area surrounding the weather station. Although satellite 
imagery can reflect the spatial and temporal distribution of precipitation, its application in various professional 
fields is limited by its poor accuracy. Fusion precipitation data is produced by fusing precipitation data from 
different sources, such as station precipitation, satellite precipitation, and radar precipitation, all on the same 
spatial and temporal scales. The advantage of combining multiple sources is that it produces data that more 
closely represents the true spatial distribution of precipitation. With the high temporal resolution CMA land data 
assimilation system (CLDAS) fusion precipitation, this study uses the El model to assess rainfall erosion in 
China on different spatial and temporal scales and to determine the potential role of rainfall on soil erosion by 
precipitation, erosive rainfall, and precipitation erosion density. The article draws the following conclusions: (1) 
CLDAS rainfall erosivity is slightly lower than the rain gauge results but presents a very good regression 
relationship with rain gauge erosivity, with high correlation coefficients at different time scales, and a lower 
margin of error than the rainfall erosivity from the Climate Prediction Center Morphing (CMORPH) technique, 
which can accurately reflect the rainfall erosivity variation on a national scale. (2) In different rainfall zones, 
trends in rainfall erosivity, rainfall volume, and the number of erosive rainfall events from 2001 to 2020 are 
generally consistent, with sharp inter-annual fluctuations in high rainfall zones. (3) Spatially, Chinese rainfall erosivity 
is characterized by high values in the southeast coastal areas and low values in the northwest inland areas. 
Temporally, erosive rainfall is concentrated from May to August. Rainfall in summer and autumn results in 
greater erosive impacts on the soil. (4) The quantitative analysis of annual rainfall, annual erosion density, and 
annual storms shows that storm and rainfall erosion density are positively correlated. In other words, the more 
rainstorm events there are, the higher the rainfall erosion density and rainfall erosivity. This study can serve as a 
scientific reference for theoretical research on soil erosion and soil conservation practices in China. 
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